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Abstract. In this paper we show that incrementally stable nonlinear time— delay systems admit symbolic 
models which are approximately equivalent, in the sense of approximate bisimulation, to the original system. 
An algorithm is presented which computes the proposed symbolic models. Termination of the algorithm in a 
finite number of steps is guaranteed by a boundedness assumption on the state and input spaces of the system. 



1. Introduction 

Symbolic models have been the object of intensive study in the last few years since they provide a tool for 
mitigating complexity in the analysis and control of large scale systems jEFPOGj . In particular, they enable 
a correct-by-design approach to the synthesis of embedded control software, see e.g. |TP06l [Tab08| . The key 
idea in this approach is to regard the synthesis of software as a control problem to be solved in conjunction 
with the synthesis of the control algorithms. Central to this approach is the possibility to construct symbolic 
models that approximately describe continuous control systems. Symbolic models are abstract mathematical 
models where each symbolic state and each symbolic label represent an aggregation of continuous states and an 
aggregation of input signals in the original continuous model. Many researchers have recently faced the problem 
of identifying classes of control systems admitting symbolic models. For example, controllable linear control 
systems and incrementally stable nonlinear control systems were shown in |TP06j and respectively in }PGT08j . 
to admit symbolic models. In the work of ,BH06j symbolic models for multi-afline systems are proposed and 
benefits from their use in solving control problems arising in systems biology and robot motion planning have 
been shown in fBBWOSj and in [BIP05], respectively. Nonlinear switched systems have been considered in 
[GPTOSj and applications to the digital control of the boost DC-DC converter have been investigated. One 
challenge in this research line is to enlarge the class of systems admitting symbolic models. In this paper we 
make a further step along this direction by focusing on the class of time-delay systems. Time-delay systems 
are an important class of dynamical systems, arising in many application domains of interest ranging from 
biology, chemical, electrical, and mechanical engineering, to economics (see e.g. jNicOH lRic03| lTbS07j ). 
In this paper we generalize the results of the work in [PGTOSj from nonlinear control systems to nonlinear 
time-delay systems. The main contribution of this paper lies in showing that incrementally stable nonlinear 
time-delay systems do admit symbolic models. An algorithm is presented which computes the proposed 
symbolic models in a finite number of steps, provided that the sets of states and inputs of the time-delay 
system are bounded. In addition to theoretical relevance, the importance of this result resides in the offer of 
an alternative design methodology based on symbolic modelH to the control design of nonlinear time-delay 
systems, which is at the present a difficult task to deal with, by using current methodologies 'Ric03'. 
In the following we will use a notation which is standard within both the control and computer science 
community. However for the sake of completeness, a detailed list of the employed notation is included in the 
Appendix (Section l7.1|l . 



This work has been partially supported by the Center of Excellence for Research DEWS, University of L'Aquila, Italy and by 
the National Science Foundation CAREER award 0717188. 

'^See e.g. ITP06I ITab08l for symbolic models— based control of linear and nonlinear control systems. 
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2. Time-Delay Systems 

In this paper we consider the foUowing nonhnear time-delay system: 

( x{t) = f{xt,u{t~r)), te]R+,a.e. 
^ ' \x{t)=^o{t), te[-A,0], 

where A e Rq is the maximum involved state delay, r G is the input delay, x{t) d X C M", xt G A" C 
C°{[—A,0];X), u{t) eU C R™ is the control input at time t e [— +oo[ , G is the initial condition, / 
is a functional from X x U to X. We denote by U the class of control input signals and we suppose that U 
is a subset of the set of all measurable and locally essentially bounded functions of time from [— r, +00 [ to U. 
Moreover we suppose that / is Lipschitz on bounded sets, i.e. for every bounded set K C X x U, there exists 
a constant k > such that 

- /(X2,U2)|| < k{\\xi - X2II00 + H^i - ""211), 

for all (xi, ui), {x2, U2) & K. Without loss of generality we assume /(0, 0) = 0, thus ensuring that x{t) = is 
the trivial solution for the unforced system x{t) — f{xt,0). Multiple discrete non-commensurate as well as 
distributed delays can appear in (j2.ip . Assumptions on / ensure existence and uniqueness of the solutions of 
the differential equation in ()2.1|) . In the following x{t,^o,u) and 2^4(^0, u) will denote the solutions in X and 
respectively in X, of the time-delay system with initial condition ^0 and input u G Z//, at time t. A time-delay 
system is said to be forward complete if every solution is defined on [0, +oo[. In the further developments we 
refer to a time-delay system as in (j2.ip by means of the tuple: 

i] = (x,A',eo,c/,w,/), 

where each entity has been defined before. 

3. Incremental Stability 

The results presented in this paper will assume certain stability assumptions that we introduce in this section. 
The following definition has been obtained as a natural generalization of the one in |Ang02| . 

Definition 3.1. A time-delay system E = {X, X,£^q,U,U, f) is incrementally input-to-state stable (J-ISS) if 
it is forward complete and there exist a ICC function (3 and a K. function 7 such that for any time t G M.q , any 
initial conditions ^1,^2 G and any inputs ui,M2 G ^ the following inequality holds: 

(3.1) ||a^t(Ci,ui)-a;t(6,ii2)IL </5(||ei-6lloo'^)+7(||K-"2)|[-.,t-.)L). 

The above definition can be thought of as an incremental version of the notion of input -to-state stability (ISS) . 
In general, inequality in (|3.ip is difficult to check directly. We therefore provide hereafter a characterization 
of (5-ISS, in terms of Liapunov-Krasovskii functionals (see |PJ06| . as far as the ISS is concerned). 

Definition 3.2. Given a time-delay system S — {X,X,£^q,U,U, f ), a locally Lipschitz functional V : 
C°([-A,0];R") X C°([-A,0];M") R+ is said to be a (S-ISS Liapunov-Krasovskii functional for E if there 
exist /Ctx3 functions ai, a2 and /C functions as, p such that: 

(i) for all a;i,a;2 G CO(hA,0];M") 

ai(||a;i(0)-a;2(0)||) < VixuX2) < a2{Ma{xi - X2)), 
where Ma : C'^([— A, 0]; R") R+ is a continuous functional such that 

7^(||x(0)||) < Ma{x) < 7a(lkl|oo), Vx G C° ( [- A, 0] ; R" ) , 
for some /Coo functions 7 and 7^; 
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(ii) for all x^.x^ G C"([-A,0];R") and ui,U2 G K" for which Ma{xx - x^) > p(||ui - U2\\) the following 
inequality holds: 

D + V{xi,X2,Ui,U2) < -a3{Ma{xi ~X2)), 

where D'^V{xi,X2,ui,U2) is the derivative of functional V in the formulation proposed by Driver 
|Dri62j . i.e. 

D^V(xi,X2,ui,U2) = limsup , 

where xf{s) = Xi{s + 0), if s e [-A,-0[ and xf{s) = x,{0)+is + e)f{x^,Ui), if s G [-0,0]. 
Theorem 3.3. A time-delay system S is S-ISS if it admits a 6-ISS Liapunov-Krasovskii functional. 

Proof. As pointed out in |Ang02| , the same lines of the proof used by Sontag for ISS, used also for time-delay 
systems in |P J06) . can be used here. Briefly, by results in |Pep07b| |Pep07a| , let 0i,02 € C^([— A, 0]; R") be 
a pair of initial conditions, mi,M2 a pair of input functions. Let — M2II00 = v. It can be proved that the 
set S = {(V'1,'02) e C"([-A,0];R") x CO([-A, 0];R") : V{iJi,ij2) < a2 o p{v)} is forward invariant, i.e., if 
{xto{(f>i,ui), Xto{4>2,U2)) & S for some to G , then {xtifpijUi), Xt{(j)2,U2)) G S for all t >tQ. In the interval 
[0,fo) with to G M"*", where {xt{4>i,ui),xt{4>2,U2)), eventually, does not belong to 5, the inequality in (ii) 
holds, which results for w{t) = V{xt{4>i,ui), Xt{(j)2,U2)) in D^w{t) < —a^ o a2^{w{t)) a.e., from which, by 
inequalities in (i), the following inequality holds, for a suitable ICC function /3, 

\\x{t,(f>i,Ui) - x{t,(f)2,U2)\\ < a^^ O /3{a2 O7a(ll01 - '/'2||oo),0- 

By the result concerning the set S, the following inequality holds: 

\\x(t, - x{t, 02, "2)11 < ° /3(a2 o 7a(ll<^i " '/'2||oo),t) + a^^ o 02 o p{v). 

From the above inequality, one gets: 

\\xt{4'l,Ul) ~ Xt{(j)2,U2)\\oc < e"'^*"^_^||01 - 02||oo 

+ aj-io^(a2O7^(||0i -02||oo),max{O,i-A}) 
+ a^^ o a2 o p{v) 

and by causality arguments, the inequality in (|3.ip is proved. □ 

At the present it is not known whether existence of (5~ISS Liapunov-Krasovskii functional is also a necessary 
condition for a time-delay system to be i5-ISS. Sufficient and necessary conditions for a time-delay system to 
be ISS, in terms of existence of ISS Liapunov-Krasovskii functionals can be found in [KPJ08| . 



4. Symbolic Models and Approximate Equivalence 

In this paper we use transition systems as abstract mathematical models of time-delay systems. 
Definition 4.1. A transition system is a sixtuple T = {Q, qo, L, «- , O, H), consisting of: 

• A set of states Q; 

• An initial state qo G Q; 

• A set of labels L; 

• A transition relation C Q x L x Q; 

• An output set O; 

• An output function H . Q O. 

A transition system T is said to be: metric, if the output set O is equipped with a metric d : O x O ^ M, 
countable, if Q and L are countable sets; finite/symbolic, if Q and L are finite sets. 



4 



GIORDANO POLA, PIERDOMENICO PEPE, MARIA D. DI BENEDETTO AND PAULO TABUADA 



We will follow standard practice and denote an element {q,l,p) £ >- by q » p. Transition systems 

capture dynamics through the transition relation. For any states q,p £ Q, q — — ^ p simply means that it is 
possible to evolve from state q to state p under the action labeled by L In this paper we will show how to 
construct symbolic models that are approximately equivalent to E. The notion of equivalence that we consider 
is the one of bisimulation equivalence |Mil89[ IParSl] . Bisimulation relations are standard mechanisms to relate 
the properties of transition systems. Intuitively, a bisimulation relation between a pair of transition systems 
Ti and T2 is a relation between the corresponding sets of states explaining how a state trajectory si of Ti can 
be transformed into a state trajectory S2 of T2 and vice versa. While typical bisimulation relations require 
that Si and S2 are observationally indistinguishable, that is -ffi(si) — H2{s2), we shall relax this by requiring 
Hi{si) to simply be close to H2{s2) where closeness is measured with respect to the metric on the output set. 
The following notion has been introduced in [GP07| and in a slightly different formulation in |Tab08| . 

Definition 4.2. Let Ti = (Qi,q^, Li, — ^-^ ,0,Hi) and T2 = {Q2,q2, L2, — ^ , 0,i?2) be metric transition 

systems with the same output set O and metric d, and let e E be a given precision. A relation i? C Qi x Q2 
is said to be an s- approximate bisimulation relation between Ti and T2, if for any (gi, q2) £ R- 

(i) diHi{qi),H2{q2))<e; 

(ii) qi — ^ pi implies existence of q2 — ^ P2 such that {pi,p2) G R] 

(iii) q2 —J* P2 implies existence of qi — -j-* Pi such that {pi,p2) G R- 
Moreover Ti is said to be e~bisimilar to T2 if: 

(iv) there exists an e-approximate bisimulation relation R between Ti and T2 such that (q^, ^ 

5. Approximately Bisimilar Symbolic Models 

In this paper we consider time-delay systems with digital controllers, i.e. time-delay systems where control 
inputs are piecewise-constant. In many concrete applications controllers are implemented through digital 
devices and this motivates our interest for this class of control systems. In the following we refer to time-delay 
systems with digital controllers as digital time-delay systems. From now on we suppose that the set U of input 
values contains the origin and that it is a hyper rectangle of the form U := [ai,6i] x [a2,62] x ... x [am,^m], 
for some < fe;, i = 1, 2, m. Furthermore given r e M+, we consider the following class of control inputs: 

, , ( u £ U : the time domain of u is [— r, — r + t] 1 

' :- I and M(t) = u(-r),t e [-r,-r + r] J' 

Given k € R" we denote by Uk.r the class of control inputs obtained by the concatenation of k control inputs 
in Ur- Given a digital time-delay system S define the transition system T^(S) := (Qi, gj, Li, — ^ , Oi, Hi), 

where: 

• Li = {li G Ur I Xr{x, li) is defined for all x G X}; 

• q -7* P, if Xr{q,li) =p; 

• Hi = Ix. 

Transition system Tr{T,) can be thought of as a time discretization of S. Transition system Tr{T,) is metric 
when we regard Oi — X as being equipped with the metric d{p, q) = \\p — q\\oo- Note that transition system 
Tt-(E) is not symbolic, since the set of states Qi is a functional space. The construction of symbolic models 
for digital time-delay systems relies upon approximations of the set of reachable states and of the space of 
input signals. Given a digital time-delay system E let RtC^) C X he the set of reachable states of S at times 
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t = 0,T, kr, i.e. the collection of all states x G A" for which there exist fc e N and a control input u G U^.t 
so that X = Xfcr(CojM). The sets i?T(S) and Ur, corresponding tc0 Qi and Li in Tr{T,) are functional spaces 
and therefore are needed to be approximated, in the sense of the following definition. 

Definition 5.1. Consider a functional space y C C^(I,Y) with Y C M", / = [a,b], a, 6 G M, a < 6. A map 
A : M"*" '2p is a countable approximation of y if for any desired precision A G M^: 

(i) ^(A) is a countable set; 

(ii) for any ?/ G 3^ there exists z G ^(A) so that \y — z||oo < A; 

(iii) for any z G .4(A) there exists y G 3^ so that \y — zjloo < A. 

A countable approximation Au of IAt can be easily obtained by defining for any \u G , 

(5.2) Au{\u) - {u e : u{t) = u(-r) G \U\2Xu , t e [-r, -r + r]}, 

where [U]2\u is defined as in (|7.ip . By comparing Z^t- in (|5.ip and .4iy(Aiy) in (|5.2p it is readily seen that 
•Aui^u) C Ur for any \u G R+. Under assumptions on U, the set .4i:^(A;^) is nonemptjH for any \u G M+. 
The definition of countable approximations of the set of reachable states RtC^) is more involved since Rri^) 
is a functional space. Let us assume as a first step existence of a countable approximation Ax of i?r(S). (In 
the further development we will derive conditions ensuring existence and construction of Ax-) We now have 
all the ingredients to define a countable transition system that will approximate Tt-(S). Given any r G K"*", 
Xx G and Xu G define the following transition system: 

(5.3) T,,A;„a„(S) := iQ2,qlL2, — ,02,i?2), 
where: 

• Q2 = .4a'(Aa'); 

• f?" S Q2 so that llfo - <72lloo < Xx; 

• L2 = .4iy(A^,); 

• q P, if lb - Xr{q, OIL ^ ^x; 

• 02 = X; 

• H2 = t: Q2 ^ 02- 

Parameters Xx and Aj^ can be thought of as quantizations of the set Rt (E) and of the space Ur , respectively. 
By construction, transition system in (|5.3p is countable. We can now state the following result that relates 
(5-ISS to the existence of symbolic models for time-delay systems. 

Theorem 5.2. Consider a digital time-delay system E = {X, X, U, lAr, f) o,nd any desired precision e G M"*" . 
Suppose that E is S-ISS and choose r G so that /3(e,T) < e. Moreover suppose that there exists a countable 
approximation Ax of Rt{T,). Then, for any Xx & K"*" and Xu G M"*" satisfying the following inequality: 

(5.4) P{e,T)+-i{Xu) + Xx<e 
transition systems Tr.x^.Xi^i^) and Tr{^) are e-bisimilar. 

Proof. The proof can be given along the lines of Theorem 5.1 in jPGT08| . We include it here for the sake of 
completeness. Consider the relation R C Qi x Q2 defined by {x, q) £ R if and only if ||ffi(a;) — ^^2(9)1100 < £• 
We now show that R is an e-approximate bisimulation relation between Tt-(E) and Tt-^Aa-.a^ (S). Consider any 
{x,q) G R. Condition (i) in Definition 14.21 is satisfied by the definition of R. Let us now show that condition 



^In fact the set Qi of states of Tt-(S) is X and not i?T-(S). However, all states in A'\i?T-(E) will be never reached and this is 
the reason why we will approximate i?r{S) rather than X. 



•^For any Xu G R+ the set An{\n) contains at least the identically null input function. 
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(ii) in Definition 14.21 holds. Consider any li G Li and the transition x — ^ y in T,-(S). By definition of L2 
there exists I2 G L2 so that: 

(5.5) II/1-/2IL <Ai,. 

Set z — Xriq, h)- Note that since I2 & L2 Q Ur, function z is well defined and z G Rr{T,). By definition of Q2 
there exists p £ Q2 so that: 

(5.6) Ik-PIL<A^. 

By the above inequality it is clear that q — ^ p in Tr^x^^\i^{Y^). Since E is S-ISS and by (|5.4p . (|5.5p and (|5.6p . 
the following chain of inequalities holds: 

||y -pIIoo = \\y - z + z - p\\oo < \\y- z\\oo + \\z-p\\oo 
< P{\\x - q\\oo,T) + 7(||/i - /alloc) + Xx 

(5.7) <P{e,T)+-f{\u) + Xx<e. 

Hence {y,p) G R and condition (ii) in Definition 14.21 holds. Condition (iii) can be shown by using a similar 
reasoning. Finally by the inequality in (|5.4p and the definition of \\£,q — 92 II ^ < £ and hence, condition 
(iv) is also satisfied. □ 

The above result relies upon existence of a countable approximation for the set of reachable states. In order 
to address this issue, we consider one possible approximation scheme of functional spaces based on spline 
analysis |Sch73j . Spline based approximation schemes have been extensively used in the literature of time- 
delay systems (see e.g. IGMPOOj and the references therein). Let us consider the space y C C°(/, F) with 
Y C R", / = [a,b], a,b eR and a < b. Given N consider the following functions (see [Sch73| ): 

J 1 — (< — a)//i, tC^[a,a + h], 
\ 0, otherwise. 



so{t) 



l — i + (t~a)/h, t [a + (i — l)h, a + ih], 
i{t) = { 1 + {t-a)/h, te [a + ih,a+{i + l)h], 

0, otherwise, i = 1, 2, A^; 

r l + {t~b)/h, te[b-r,b], 
'^+1^^^ - \ 0, otherwise, 

where h = {b — a)/{N+l). Functions Si called splines, are used to approximate y. The approximation scheme 
that we use is composed of two steps: we first approximate a function y G 3^ (Figure [TJ upper panel) by means 
of the piecewise-linear function yi (Figure [U medium panel), obtained by the linear combination of the + 2 
splines Si, centered at time t = a + ih with amplitudeo y{a + ih); we then approximate function yi by means 
of function 7/2 (Figure [TJ lower panel), obtained by the linear combination of the A'^ + 2 splines s^, centered 
at time t = a + ih with amplitude jji in the latticql [^lae, which minimizes the distance fron{^ y(a + ih), i.e. 
iji = argminj^g[y]2^ ||y - y{a + ih)\\. 

Given any G N, 61, M G K+ lelQ: 

(5.8) A{N,e,M) := h'^M/8 + {N + 2)9, 

with h ^ [b — a)/iN + 1). Function A will be shown to be an upper bound to the error associated with the 
approximation scheme that we propose. It is readily seen that for any A G and any M G there always 



^This first step allows us to approximate the infinite dimensional space y by means of the finite dimensional space V^"*"^. 
^We recall that the set [Y]2e is defined as in 117.11 1. 

''This second step allows us to approximate the finite dimensional space jjy means of the countable space ([i^]2e)^^'^i 

which becomes a finite set when the set Y is bounded. 

^The real M is a parameter associated with y and its role will become clear in the subsequent developments. 
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Figure 1. Spline-based approximation scheme of a functional space. 



exist N (z N and 9 G M^" so that A(N,9, M) < A. Let A^a,a/ Eind 6\^m be such that A{Nx^Mt^\,Mj 

M) < A. 

For any A G K+ and M G IR.+ , define the operator 

iJx.M ■y^C°{[a,b];Y), 
that associates to any function y £y the function: 

(5.9) i'x.MivW) ■■= msr{t), ie[a,6], 

where yi G [i^]2eA j\f ^-^^^ \\yi — y{a + < 9\^m, for any i = 0, 1, Nx^m + 1- Note that operator ipx.M is 
not uniquely defined. For any given M G M"*" and any given precision A G define: 

(5.10) Ay,MW-=^KM{y). 

The above approximation scheme is employed to construct countable approximations of the set _R,-(E) of 
reachable states (see Proposition 15. 3p . 

Consider a digital time-delay system E = (X, X , £^Q,U,Ur, f) and suppose that: 

(A.l) Eis5-ISS; 

(A. 2) X and U are bounded sets; 

(A.3) Functional / is Frechet differentiable in C°([-A, 0]; M") x R"'; 

(A.4) The Frechet differential J(0, u) of / is bounded on bounded subsets of C°([-A, 0]; R") x R™. 
Under the above assumptions, the following bounds are well defined: 

Bx = snp^^x W^t Bj = sup(^,„)gcn([_A,o];X)x;7 \\Ji(t>,u)\\, 
Bu ^sup^^uM, M = {PiBx,0)+iiBu) + Bu)^Bj, 

where k is the Lipschitz constant of functional / in the bounded set C°([— A, 0]; X) x U and || J((/), m)|| denotes 
the norm of the operator J{(p,u) : C°([— A,0]; R") x R™ R". We can now give the following result that 
points out sufficient conditions for the existence of countable approximations of Rr{'S). 

Proposition 5.3. Consider a digital time-delay system S — {X, <Y, ^q, U,lAr, /), satisfying assumptions (A.l- 
4 ) and the following conditions: 



(A5) ^oePCH[-A,0];X), HCoIU < < S^, < M, 

/3(i?o.,0)+7(Sc/) <i?x, /3(i?°:,T)+7(Bc;) <i?o., r > 2A, 



8 GIORDANO POLA, PIERDOMENICO PEPE, MARIA D. DI BENEDETTO AND PAULO TABUADA 

with M as in 115. Then the set Ax defined for any Xx £ by: 

(5.11) Ax{Xx) = iJX^,M{Rr{^)), 

with 'ip\x,M CIS in I15.9\) . is a countable approximation of Rr{T,). 
input: 

time-delay system E {X, X, ^o, U,U, /) satisfying assumptions (A. 1-5); 
parameters r, iV, 6, \u,M; 
init: 
k := 0; 

Q'^ ■= {92}: where = ipx.Aii^o), with t/ja.m defined as in ([SJ]) and A = A{N,9,M); 
g'=-i 0; 



k 

H2 := i:Q2^ O2; 
h:=A/{N + l); 
while ^ g'=-i do 
foreach q E Q'^ do 

foreach I2 G [C^]2Am do 
compute z := Xr{q, I2)', 

compute p = tp\ m{z), with -0a m defined as in (|5.9p and A = A(7V, 9, M); Q'^+^ := Q'^ U {p}; 

— ^= — u{(g,/2,p)}; 

end 
end 

k:=k+l; 
end 

output: T,,jv,e,A„(S) (Q^ 9^, [C/]a„ , — ,-^,^^2) 

Algorithm 1: Construction of symbohc models for time-delay systems. 



While the first and third assumptions in (A. 5) are given apriori on the time-delay system under study, the 
other assumptions are satisfied for sufficiently large values of Bx and r. Note that if Bx does not satisfy the 
assumptions in (A. 5) one can always embed the state space X of the time-delay system in a bigger state space 
X', so that the corresponding bound Bx' = sup^^x' W^W satisfies the required assumptions. The proof of the 
above result requires some technicalities and is therefore reported in the Appendix (Section l7.2p . 
We now have all the ingredients to define a symbolic model for digital time-delay systems. Given r e R+, 
0, \u G R"*" and S N, consider the transition system 

(5.12) Tr^Nfi^ui^) (g2,'Z2°,^2, — ,02,i?2), 

where: 

• Q2 = ^A'(A(7V, e, M)) with Ax as in ([5lT|) with Xx A(Af, 6*, M) and M as in ([5lT|) : 

• ^2 = V'a,m(Co), with i'\;^.,M defined as in (|5.9p and Xx = A{N,0, M); 

• L2^ Au{Xu); 

, q^p,ii \\p-Xr{q,l)\\^ < A{N,e,M); 

• 02 = X- 

• H2=i:Q2^02. 

Note that the transition system in (|5.12|) coincides with the one in (|5.3p by setting Xx = A(A^, 9, M). Moreover, 
it is readily seen that: 

Proposition 5.4. If the digital time-delay system S satisfies assumptions (A. 1-5), transition system Tr^N,e,\u (^) 
in Ii5.12\) is symbolic. 
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Transition system TTj^ g xuC^) can be constructed by analytical and/or numerical integration of the solutions 
of the time-delay system. One possible construction scheme is illustrated in Algorithm [1] which proceeds, as 
follows. The set Q'' of states of the symbolic model at step fc = is initialized to contain the (only) symbol 
12 = V'a,m(^o) that is associated with the initial condition ^o- Then, for any initial condition q G Q'' and any 
control input I2 € [C^]2Au, the algorithm computes the solution z = Xr{q, h) of the differential equation in (j2.ip 
at time t = t, and it adds the symbol p = ipXM^^) to ■ In the end of this basic step, index k is increased to 
A: + 1 and the above basic step is repeated. The algorithm continues by adding symbols to since no more 
symbols are found, or equivalently, since a step k* is found, for which Q'^ — ^'^ ■ Convergence properties 
of Algorithm [T] are discussed in the following result. 

Theorem 5.5. Algorithm[^ terminates in a finite number of steps. 

Proof. Let Z be the collection of all functions of the form (|5.9|) with yi,2/2, •■•,2/a'a ^ i^he- Since the set 
X is bounded, the set [X]2B is finite and hence the set Z is finite as well. By construction the sequence Q'^ is 
non-decreasing, i.e. Q'' C (5'^+^ and each set of the sequence is contained in Z, i.e. Q'' C Z. Hence, a fixed 
point of Algorithm [1] will be found in a finite number of steps, which is upper bounded by the cardinality of 
Z. □ 

We can now give the main result of this paper. 

Theorem 5.6. Consider a digital time-delay system S — (X, X , ^0, U,lAr, f) and any desired precision e G M^. 
Suppose that assumptions (A. 1-5) are satisfied. Moreover let t,9,Xi( G and N gN satisfy the following 
inequality 

(5.13) /3(£,T)+7(A;^)+A(iV,^,Af) <e, 

with A as in i5.8\) and M as in i5.11\) . Then transition systems 7V(I]) and Tr^N,e,Xui^) ^'"6 e-bisimilar. 

Proof. The map Au is a countable approximation of U and by Proposition 15.31 the map Ax is a countable 
approximation of i?r(S). Choose Xx G and Xu G satisfying the inequality in (|5.4p . There exist 9 G M"*" 
and G N so that Xx = A{N,6,M) and hence the inequality in (|5.13[) holds. Finally the result holds as a 
direct application of Theorem 15.21 □ 

The above result is important because it provides a method to translate time-delay systems to approximately 
bisimilar symbolic models. Hence, it gives a concrete alternative methodology to the control design of nonlinear 
time-delay systems, by regarding control design on time-delay systems as control design on symbolic models 
(see e.g. }TP06| [Tab08] for symbolic models-based control of linear and nonlinear control systems). 

6. Conclusion 

In this paper we showed that incrementally input-to-state stable digital time-delay systems admit symbolic 
models that are approximately bisimilar to the original system, with a precision that can be rendered as small 
as desired. An algorithm has been presented which computes the proposed symbolic models. Termination of 
the algorithm in finite time is ensured under a boundedness assumption on the sets of states and inputs of the 
system. 
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7. Appendix 

7.1. Notation. The symbols N, Z, M, M+ and M.^ denote the sets of natural, integer, real, positive and 

nonnegative real numbers, respectively. Given a vector a; G the i-th element of x is denoted by Xi; 

furthermore ||a::|| denotes the infinity norm of .x; we recall that := maa;{|a;i|, \x2\, \xn\}, where \xi\ is the 
absolute value of Xi . For any A C M" and 9 G define 

(7.1) [A]0:={aeA\a^=kie, h £ Z,i = 1, ...,n}. 

Given a measurable and locally essentally bounded function / : Rq M", the (essential) supremum norm of 

/ is denoted by ||/||oc; we recall that ||/||oo := (ess)sup{||/(t)||, t > 0}. For a given time r e M+, define 
so that frit) = fit), for any t £ [0, r[, and f{t) = elsewhere; / is said to be locally essentially bounded if 
for any r G M+, fr is essentially bounded. A continuous function 7 : R'^ R'^ is said to belong to class IC 
if it is strictly increasing and 7(0) = 0; 7 is said to belong to class /Coo if 7 G /C and "f{r) 00 as r —>■ 00. 
A continuous function /3 : Rq x Rq — » Rq is said to belong to class IC£ if for each fixed s, the map /?(r, s) 
belongs to class JC with respect to r and, for each fixed r, the map /3(r, s) is decreasing with respect to s and 
l3{r,s) as s ^ 00. Given k,n e N with n > 1 and / = [a, 6] C R, a, 6 G R, a < 5 let C'=(/;R") be the 
space of functions / : / — > R" that are continuously difi'erentiable k times. Given A: > 1, let PC'^(/;R") be 
the space of C'^~^{I; R") functions / : 7 ^ R" whose fc-th derivative exists except in a finite number of reals, 
and it is bounded, i.e. there exist 70,71, •.•,7s G R"*" with a = 70 < 71 < ... < 7s = & so that D*^ / is defined 
on each open interval (7^, 7i+i), i = 0, 1, s — 1 and maxj=o,i,...,s-i ^^Ptei'yi,'yi+i) W^'^ /WIloo < 00. For any 
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continuous function x{s), defined on —A < s < a, a > 0, and any fixed t, < t < a, the standard symbol Xt 
will denote the element of C°([-A, 0]; M") defined hy Xti0) = x(t + 6), -A<0 <0. The identity map on a 
set A is denoted by 1a- Given two sets A and B, if A is a subset of B we denote hy ia ■ A ^ B ot simply by 
I the natural inclusion map taking any a ^ A to i{a) = a ^ B . Given a function f : A ^ B the symbol f{A) 
denotes the image of A through /, i.e. f{A) {6 G S : 3a G A s.t. b = /(a)}- 



7.2. Technical Proofs. The proof of Proposition 15 . 31 is based on the following lemmas. 

Lemma 7.1. Suppose that y C PC^{I;Y) and there exists M G IR+ so that WD"^ y\\ao < M for any y &y. 
Then Ay.m as defined in 115. 10\) . is a countable approximation ofy. 

Proof. Let Z be the collection of all functions of the form (|5.9p with yi,y2, ■■■,yNx m+i ^ [^126^ Af • 'con- 
struction, since i/'a,m(3^) is a subset of Z that is countable, it is countable as well. Hence, condition (i) in 
Definition 15.11 is satisfied. Let us now show that also condition (ii) is satisfied. Consider any A G and any 
y G y and set Hx^m ^ {b — cl)/{Nx^m + !)• By Theorem 2.6 in |Sch73j and the definition of M, the following 
inequality holds: 

(7.2) \\y-7r\\oo < hl,,\\D'y\U8 < hl,,M/8, 

where Tr{t) - Y.fX'^^ V (^^a.m* + «) s,{t), t G [a, b]. Moreover by setting z{t) = V'A,A/(2/(t)) = E^'o"^' V^s^it)^ 
t G [a, 5], the following chain of inequalities holds: 

< 



YIi=o '^^ [y {h\,Mi + a)-yi)Sr 



OO 



E^r^'ll(2/(W + a)-y.)s.lL< 



T^X'''\\y{h^.M^ + a)-UM\oo< 
(maxi=o,i,...,AfA M+i \\y {hx,Mi + a) - yj) J2i=o^' PJoo < 

By combining inequalities in (|7.2p and in (|7.3p and by definition of 6x^m and Nx.m, one gets: 

l|y-z||oo < lly - ttIIoo + Ik - ziloo 

< hl„iM/8 + 9x,MiNx,M + 2) - A{Nx,M, Oxm^M) < A. 

Hence, condition (ii) in Definition 15.11 is satisfied. We conclude by showing that also condition (iii) holds. 
Consider any z G Ay,M{X). By construction there exists y G 3^ so that z — '4'x,M{y)- Hence, by following the 
same reasoning in proving condition (ii), condition (iii) can be proved as well. □ 

Under assumptions in (A. 1-4), the regularity properties of the initial state in (A. 5) propagate to the whole 
set of reachable states, or in other words, time-delay systems are invariant with respect to those properties in 
(A. 5). More precisely: 

Lemma 7.2. Consider a digital time-delay system S — {X, X,£^q,U,U, f), satisfying assumptions (A. 1-5). 
Then for any Xr G i?r(S), 



x^ ePC^{[-A,0]-X), \\x^\U<B"x, ||2)VL<M. 

Proof. First note that the function t —> x{t), t G [0, t], is uniformly continuous in the (compact) set [0, r]. Since 
T > 2A, it follows that Xr+e e C^{[-A,0];X), 9 G (-A,0) (i.e. the derivative i^+e belongs to C°([-A, 0]; AT)). 
Moreover, by taking into account the Lipschitz property of /, the (5-ISS inequality, the bounds on initial state 
and input, the following inequality holds: 

||ir+e||oo = sup^^[_^ o^\\f{xr+e+a,u{T + 9 + a-r))\\ 

< KSup„g[_A,o](l|2;r+e+Q||co + \\u{t + 9 + a - r)||) 

< K{(3iB^x,0)+^{Bu) + Bu),ee]-A,0[. 
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As far as the second derivative is concerned, the foUowing equahty holds, for 6* g] — A, 0[, 



^g|2 - JiXr+0,u{T + e -r)) . Q 

By taking into accounts the bound on the Frechet differential, and the bound on the derivative Xr+e and 
u{t) — 0, we obtain H-D^Xt-Hoo < M. Finally by assumptions of -B^, Bu and t in (A. 5), it is readily seen that 
||a:.||oo <5!^. □ 



By combining Lemmas 17.11 and 17. 2i the proof of Proposition 15.31 holds as a direct consequence. 
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